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Abstract 

Recent results on the semiclassical dynamics of an electron in a solid are explained using 
techniques developed for "exotic" Galilean dynamics. The system is indeed Hamiltonian and 
Liouville's theorem holds for the symplectic volume form. Suitably defined quantities satisfy 
hydrodynamic equations. 
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The effective semiciassical dynamics of a Bloch electron in a solid is modified by a Berry 
curvature term J^. Putting h = 1, the equations of motion in the n^^ band read indeed 

r =^-kxn(k), (1) 

k = -eE - er X B(r), (2) 

where r = (x*) and k = (fcj) denote the electron's intracell position and quasimomentum, re- 
spectively, en(k) is the band energy. The purely momentum-dependent ft is the Berry curvature 
of the electronic Bloch states, ^^^(k) = eijidk_.Ai(k), where A is the Berry connection. The elec- 
tronic charge is — e. The new Berry-curvature dependent term in eqn. has been instrumental 
in explaining the anomalous Hall effect in ferromagnetic materials [2], and also the spin Hall 
effect 0. 

Much of the transport phenomena rely on the Liouville's theorem on the conservation of the 
phase-space volume element. Recently, Xiao et al. observed that, due to the Berry curvature 
term, the naive phase-space volume form drdk is no more conserved. They attribute this fact to 
the non-Hamiltonian character of the equations ©-(El), and suggest to redefine the phase-space 
density by including a pre- factor, 

Pn = Dnfn{r,Kt), D = (1 + eB • fi) , (3) 

where d is the spatial dimension and /„ is the occupation number of the state labeled by r, k. 
If the latter is conserved along the motion, dfn/dt = 0, then /)„ satisfies the continuity equation 
on phase space. Eqn. @ is the starting point of Xiao et al. to derive interesting applications. 

Xiao et al. derive Eq. 0, their main result, by calculating the change of drdk using 
the equations of motion and then finding a compensating factor. This may appear somewhat 
mysterious at first sight. Our Note is to the end to point out that the equations ©-(El are 
indeed Hamiltonian in a by-now standard sense [HI, and the validity of Liouville's theorem is 
restored if the symplectic volume form, Eq. (|7|) below, is used. Then Eq. ^ follows at once. 
Our clue is the relation to "exotic" Galilean dynamics, introduced independently and equivalent 
to non-commutative mechanics in (51 [3 |H] . 

Let us first summarize the general principles of Hamiltonian dynamics . The first ingre- 
dient is the symplectic form on 2ci-dimensional phase space, i.e., a closed and regular 2-form 
Lo = \uJapd$^'^ A d^^, where denotes the phase space coordinates r and k collectively. Then 
the Poisson bracket of two functions / and g on phase space is 



{f,g}=uj"^do.fdpg, (4) 

where a;"'^ is the inverse of the symplectic matrix, uj'^'^LO^p = (5"^. Since lo is closed, the Poisson 
bracket (HJ satisfies the Jacobi identity. Then Hamilton's equations read 

r = {/i,n, (5) 

where h = h{^) is some given Hamiltonian, which is the second ingredient we need. Such a 
framework is obtained ^1, in particular, starting with a first-order Lagrangian L = aa(^)^" — 
h{^) on phase-space, whose Euler-Lagrange equations read 

^api'^ = dah, where uJap = d^ap - dpa^, (6) 

and multiplying Q by the inverse matrix yields Hamilton's equations, ©. 



Now the natural volume form on phase space is the d*^ power of the symplectic form [Sj 
which, in arbitrary coordinates on phase space, reads 

2d 

dV = -u:^ = ^det(w„^) n (7) 

0=1 

The intrinsic expression Q is invariant w.r.t. coordinate transformations. The naive expres- 
sion n^?") changes in fact as H ~^ 11 '^C" = det(9"C/t^^^) 11 under a transformation 
^ — > C. This is, however, compensated by the change of the determinant. As a bonus, the 
validity of Liouville's theorem is restored: the r.h.s. of © generates the classical flow of the 
phase space, w.r.t. which the symplectic volume form, is invariant, since the classical flow 
is made of symplectic transformations [S] . 

Now Darboux's theorem tells that one can always find canonical local coordinates Pi^q' 
(say) such that the symplectic matrix has the canonical form to = dpi A dq^. Then the general 
Hamilton equations, © take their familiar form qi = dp.h, pj = —dq.h. The volume form is 
simply Y{dpif\dq^ . Xiao et al. use the restricted terminology, referring to canonical coordinates. 
Non-canonical coordinates may be useful, though. The equations of motion Q-© can indeed 
be put into the form © with symplectic and resp. Poisson matrix 

(—EijkeB^ —Sjk \ oil ^ijk^^ 5jk + eBjQk 

6,k ; 1 + eB • n \^ _ ^^^^k -SijueB'^ 

where i, j, A; = 1, 2, 3. The Poisson bracket relations are therefore 



{x\x^] 
{x , kj} 



' 1 eB • n ' 



The Berry curvature is divergence- free, d\<^.Vt^ = 0, because it is a curl. The Jacobi identity 
follows hence from d^jB^ = 0, cf. The Hamiltonian is h = Cn — eV. Note, in particular, that 
the coordinates do not commute: the model is indeed equivalent to non- commutative mechanics 

[sua El- 

When both the B, and the Berry curvature, ft =const. are directed in the third direction, 
these relations reduce to those of the "exotic" model in the plane, derived from first principles 
[HI IS]. The crucial anomalous velocity term in arises, in particular, by minimal coupling of 
the "exotic" model to a gauge field Eq. Q is obtained, furthermore, by the Faddeev-Jackiw 
jlUj construction mentioned above, if we start with the Lagrangian ^ E] 

^Bioch ^ _ eA,ir,t))x' - (e„(k) - eVir,t))+A\k)ki, (12) 

where V and Ai are the scalar resp. vector potential for the electromagnetic field and A is the 
potential for the Berry curvature. The conserved volume form is j2j 



/det(a;„/3) H A dx' = (l + eB fl)Y[ dh A dx\ (13) 

i 

Consistently with Darboux's theorem, canonical coordinates can be found. 



In the planar case, for example, if B is constant such that 1 + eBQ ^ 0, 

1 



Pi = Vl + eB^ ki + -eBsijQ^ 

(14) 



= X* - 4^ f 1 - Vl + eBn) e 
eB V / 



are canonical [7j, so that the Poisson bracket reads simply {F, G} = Qq.F dp.G — dQ^Gdp.F. 
(The disadvantage is that the Hamiltonian becomes rather complicated). When expressed in 
terms of canonical coordinates, the four-dimensional volume form reads simply dQi A dQ2 A 
dPif\dP2\ the 2D version of is recovered as dV = det(9(P, Q)/(9(k, r)) drdk. Similar formulae 
would work in higher dimensions. 

Let us consider a particle distribution / = /(r,k, t). The probability of finding our particle 
in the phase space volume element dV is, see (|13j) . 



fdV = f^dei{ujap) n di"' = /(r, k, t)(l + eB • Q,)dxidx2dx^dkidk2dk'i. (15) 

a 

Since dV = a;^/6 is invariant by the Liouville theorem, fdV is an integral invariant of the flow, 
provided df /dt = 0, i. e., dtf + {/, h} = 0. However, 

SO that / can not be considered as a density. Multiplication with the square root of the deter- 
minant removes the unwanted term, though : 



P = /y'det(wa/3) (16) 



does satisfy the continuity equation dtp + • (Cp) = on phase space (and time). In the 
coordinates (r,k), ()16() reduces to the expression ^ of Xiao et al., which is hence simply the 
coefficient of 11 (|15)) . written in a local coordinate system. The (modified) Boltzmann 

transport equation for p is 

^ = {p,h} + dtp= —^^—\(Vk€n + eExn + eB{n-Vk€n)) -VrP 
dt 1 + eB • i i L V / 

-e(E + Vke„ X B - en(E • B)) • V^p] = (17) 

where (B • Q)' = dt(B ■ Q) + Q^r • Vr-Bj + -Bjk • Vk^^* is the material derivative. 

Using the equations of motion and Maxwell equations V • B = and dtB + Vr x E = 0, a 
tedious but straightforward calculation shows, furthermore, that 

J pdk and v = — J rpdk (18) 

[where the integration is over the first Brillouin zone] satisfy the continuity equation dt^p + Vr • 
{vif) = in ordinary space-time. Similarly, we also have the Euler equation 

^{dt + V • Vr)v = f - VCT, (19) 



where f = Jrpdk is the mean force and a = (aij) with cjy = J{xi — Vi){xj — Vj)pdk is the kinetic 
stress tensor jj]. 



All these remarks hold when 

det{ujap) = (1 + eB • flf / 0. (20) 

In the singular planar case eBO, = — 1 [which can only happen if and B are constants], the 
volume element (jTj) degenerates and the phase space looses 2 dimensions. The symplectic matrix 
can not be inverted and one can not get from a Lagrangian to the Hamiltonian framework. The 
only allowed motion follow the Hall law [HI 13 El- Hamiltonian reduction yields, however, a 
well-behaved 2-dimensional phase space with symplectic form proportional to the surface form 
in the plane, LOred = —eBdx A dy. The fluid moves collectively, governed by a generalized Hall 
law 0. 
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